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Abstract 

We show that a recently developed method for generating bounds for the 
discrete energy states of the non-hermitian —ix^ potential (Handy 2001) is 
applicable to complex rotated versions of the Hamiltonian. This has im- 
portant implications for extension of the method in the analysis of resonant 
states, Regge poles, and general bound states in the complex plane (Bender 
and Boettcher (1998)). 
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I. INTRODUCTION 



In a recent work, Handy (2001) presented a novel quantization formalism for generating 
converging bounds to the (complex) discrete spectra of non-hermitian, one dimensional, 
potentials. The first part of this analysis makes use of the fact that the modulus squared of 
the wavef unction, S{x) = |^'(a;)|^, for x e 3?, and any (complex) potential function, V{x), 
satisfies a fourth order, linear differential equation: 




where S^^'^x) = d^S{x), and E^Er + iEj, etc. 

The bounded {L'^) and nonnegative solutions of this differential equation uniquely cor- 
respond to the physical states. For the case of rational fraction potentials, there is an 
associated (recursive) moment equation, whose solutions are parameterized by the energy, 
E. Constraining these through the appropriate Moment Problem conditions for nonnega- 
tive functions (Shohat and Tamarkin (1963)), generates rapidly converging bounds on the 
discrete state energies. This entire procedure is referred to as the Eigenvalue Moment 
Method (EMM) and was originally developed by Handy and Bessis (1985) and Handy et 
al (1988a,b). An important theoretical and algorithmic component is the use of linear pro- 
gramming (Chvatal (1983)). 

The EMM procedure is not solely a numerical implementation of some very important 
mathematical theorems. It is possible to implement it algebraically, deriving bounding 
formulas for the eigenenergies (Handy and Bessis (1985)). Furthermore, it does not depend 
on the hermitian or non-hermitian structure of the Hamiltonian. What is important is that 
the desired solutions be the unique nonnegative and bounded solutions of the differential 
system being investigated. In this sense, it provides an important example of how positivity 
can be used as a quantization procedure. 

Many important problems can be analyzed through the application of EMM on Eq.(l), 
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and its generalization as outlined in Sec. II. This includes generating bounds to Regge poles 
(which first motivated Handy's formalism), to be presented in a forthcoming work (Handy 
and Msezane (2001)). Other problems include the calculation of resonant states as well 
as bound states in the complex plane as discussed by Bender and Boettcher (1998), and 
elaborated upon by others cited in the more recent work of Mezincescu (2001). 

However, the realization of these depends on understanding the complex rotation exten- 
sion of Handy's original formulation. Indeed, such concerns introduce a new class of EMM 
problems not encountered before, and thereby motivate the present work. 

Our basic objective is to extend the 5'(a;)-EMM formulation to complex rotated trans- 
formations of the Hamiltonian — ix^. We are particularly interested in the effectiveness of 
such an analysis. Indeed, once the correct moment equation is developed, one finds that the 
EMM approach works very well in this case, and yields the anticipated result, as detailed in 
Sec. III. 

In the following sections we present a more general derivation of the S{x) equation, and 
apply the resulting formalism to the rotated —ix^ Hamiltonian. 
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II. GENERALIZED 5-EQUATION 



Consider the (normalized) Schrodinger equation 

- dl^(x) + V{x)^(x) = E^(x), (2) 

for complex energy, E, and complex potential, V{x). Assume that the (complex) bound 
state, ^(x), lies within the complex-x plane, along some infinite contour, C. Let x{^) define 
a difi^erentiable map from a subset of the real axis to the entire complex contour: 

x{0 -.ie^^C. (3) 

The transformed Schrodinger equation is 

- (l>(09^)'*(0 + ViO^iO = E^iO, (4) 
where D{^) = {d^x)~^, and V{^) = V{x{^)). Alternatively, we may rewrite the above as 

H^^io = mdim + md^'^io + cmio = o, (5) 

where A(^) = -(L'(0)^ ^(0 = -i5?(^(0)^ and C(0 = ^(0 - E. In the previous work 
by Handy (2001), he showed how such an expression, without the term, leads to a fourth 
order equation for S{^) = |^'({)p. We broaden this derivation to include the above general 
case. 

Define the following four functions: 

Ei(O = **(O^«*(O + c.c. = 0, (6) 

E2(6^^'*(6i/«*(0 + c.c. = 0, (7) 

Ai(O = **(e)i^e*(O-c.c. = 0. (8) 

and 

A2(0 = *'*(0^^€*(0-c.c. = 0. (9) 



Define 



s{0 = **(0*(0 



(10) 



p{0 = 



(11) 



^(0 



2i 



(12) 



and 



a^^(0a|^*(0-ag^*(09|^(0 

2i 



(13) 



These correspond to important physical quantities. The first two are nonnegative functions 
corresponding to the probabihty density and the "momentum density", while J(^) is the 
probability flux. 

We then have (i.e. A — An + iAj, etc.): 



(14) 

^ = {S"{0 - 2P(e))^7(0 + S'iOBiiO + 25(0^^/(0 - ^{Br{0 + Anm))J{0 = 0, 

(15) 

E2(0 = P'iOMO + mOMO + mOMO + s'(OCr(0 - 2 J(OC/(0 = 0, (16) 



= p'iOMO - moMO + mmio + s'(oci(o + 2 j(oc«(o = 0. (i?) 



For the case of the ordinary (real potential) Schrodinger equation where x — ^, and 

Aj = 0,B = 0, and Cj = 0, we can combine Eq.(14) (divided by A^) and Eq.(16) to 
generate a third order differential equation for 5" (Handy (1987a,b), Handy et al (1988c)). 



and 
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If A/ = and S = 0, the case considered by Handy (2001), then Eq.(14) (divided by 
Ar) yields an equation for P. Differentiating this relation, and inserting Eq. (16) defines a 
third order equation in S, which also depends on J. We can elliminate J by differentiating 
(after dividing out ^) one more time, and substituting Eq.(15). The resulting expression 
corresponds to that in Eq.(l), assuming Ar — 1. 

U Aj = and C/ ^ 0, or Ar = and Cj ^ 0, we can determine J from Eq.(16) or 
Eq.(17), respectively, and procede with the analysis discussed below. 

If both Arj 7^ 0, multiplying Eq.(16) and Eq.(17) by Aji(^) and Aj{^), respectively, and 
adding, yields an equation for J: 

-2Im{A{0C*{0) • ^ ' 

We can now substitute this expression, J(^; P, P', 5"), into Eq.(14) and Eq.(15) yielding 
two coupled differential equations for P and S. Instead of the formalism presented by Handy 
(2001), we can work with these two coupled equations, since they involve two nonnegative 
configurations (for the physical solutions). However, depending on the nature of the function 
coefficients (i.e. A{^), B{^), C{^)) it may be expedient to reduce these to one fourth order, 
linear differential equation for S. This is particularly the case corresponding to that discussed 
by Handy (2001). 

For completeness, we note that the analytic extension of S{^) onto the complex-^ plane 
(which is of no immediate concern to us) corresponds to 

^(0 = **(r)*(0- (19) 

That is, if we assume that ^(^) is analytic at the origin (for simplicity), then ^(^) = 
J^jLoCj^^ and \l/*(^*) = J2jLoC*C'' are convergent power series, defining an analytic product 
which is nonnegative along the real ^-axis. 



6 



III. THE COMPLEX ROTATED + {IX)^ HAMILTONIAN 



Consider the non-hermitian potential problem 

- ^"(x) - ix^^{x) = E^f(x), (20) 

first considered in the work by Bender and Boettcher (1998). Prom their analysis, it is 
known that bound states exist along the real axis. They satisfy 

JdxPjx) - iJdxx^Sjx) 
Jd^) • 

If the physical solution conforms with the VT invariant nature of the Hamiltonian, 

^*{-x) = ^(x), (22) 

then S{—x) = S{x), and the discrete state energy is real. Exhaustive investigations by many 
individuals supports this. These include the works by (in addition to those cited) Bender 
et al (1999), Bender et al (2000), Bender and Wang (2001), Caliceti (2000), Delabaere and 
Pham (1999), Delabaere and Trinh (2000), Levai and Znojil (2000), Mezincescu (2000,2001), 
Shin (2000), and Znojil (2000). The analysis of Handy(2001) has also been extended to the 
case of complex E, also supporting that the discrete states of the —ix^ potential do not 
violate VT symmetry (Handy (2001)). In accordance with all of this, we procede assuming 
that the discrete spectrum is real. 

Let X = e^^\x\. Adopting Bender and Boettcher 's notation, it is known that the bound 
states decrease exponentially fast in any asymptotic direction within the 'wedges' defined 

by 

\o-eL,R\<l, (23) 

where 9l = — tt + y^, and 9r = The anti-Stokes lines (along which the asymptotic 

decrease is fastest) is centered at 9l,r- 

We first considered the general transformation 
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x{0 = e''C 



e 



(24) 




for < ^ < cxo, and lying within the boundedness wedge(s). Such transformations would 
be relevant not only for the —ix^ potential, but all the other non-Hermitian potentials con- 
sidered by Bender and Boettcher. However, the generation of a moment equation, following 
the formalism in Section II, made this approach inefficient, due to the large order of the 
associated moment equation relation. 

Instead of this approach, the —ix^ problem affords a particularly simpler formulation, 
since the wedges intersect the entire real axis. This is not the case for all the other potentials 
(i.e. (— ix)", n > 4) considered by them. The generalization of the present approach can 
be made to these cases as well; however, we defer this to another communication, since it 
requires a completely different linear programming formalism. 

Consistent with the formalism in Sec. II, if we take 



where ^ e (— oo,-|-oo), then the bound state solution ^(e*^^) will be asymptotically zero 
(i.e. at ^ = ±oo) for 

-ro<'><ro- 

Clearly, this includes the real axis. However, for fixed 0, the closer e*^^ is to the anti- 
Stokes line along the positive real-^ axis (i.e. 9 Or) the farther is e*^^ from its respective 
anti-Stokes line, for ^ < (and conversely). Only along the real axis {9 = 0), will the 
exponential decay be equally balanced between both the positive and negative real axis. 

The work by Handy(2001) introduces a very accurate bounding technique which was 
used to generate the first five low lying states. Our immediate interest is in testing the same 
EMM formalism with respect to complex rotations of the system. 




(25) 
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A. Moment Equation for the Complex Rotated —ix^ Potential 



As noted, we will be solving the Schrodinger equation along the complex- a; infinite con- 
tour defined hy x = e*^^, for — oo < C < +00. The resulting differential equation is 



- - ^e''^^'*(0 = Ee^^^-^{^). (27) 

We will assume E to be real. The resulting, fourth order differential equation for S{^) is 
then 

where 

A{0^c,e + Es2, (29) 

Cn = cos{ne), Sn = sin{ne) , and ^("^(O = d^S{0. 

We know that all of the solutions to this equation are regular in the entire complex ^ 
plane, despite the singular function coefficient, (Handy(2001)). There is a real singularity 
at ^ = 6, defined by 

6=-(±^)3eSR. (30) 
C5 

In terms of the implicit 9 dependence, lim^^o K^) = 0. 

In all of the previous applications of the EMM approach, we have never considered prob- 
lems with function coefficients singular within the physical domain. The types of singularities 
encountered have either been close to the physical domain, or at the boundaries. 

Despite the fact that the singularities in question do not affect the regular nature of the 
solutions, they are of concern within the EMM formalism, as discussed in the recent work by 
Handy, Trallero-Giner, and Rodriguez (HTR, 2001). Intuitively, let us represent the above 
fourth order equation by O^S{^). If one is to generate a moment equation by multiplying 
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both sides of this differential equation by A(^)^, one must make sure that the resulting 
moment equation does not admit solutions to the more general problem of O^S{^) = V{^—b), 
where — b) is a distribution- like expression supported at b. Such expressions disappear 
when multiplied by A(^) (with a zero at b). If one is not careful, the naive moment equation 
may not correspond to the desired problem, V — 0. That is, an inappropriate specification of 
the moment equation would not produce any converging bounds. Such concerns are generally 
inconsequential, within the EMM framework, if it is known that the desired solution also 
has a zero at b (HTR, 2001). For the present problem, this is not the case, and a suitable 
modification of the standard EMM formalism is required. 

In order to address this problem, one should first perform a translation change of vari- 
ables: 

X-^-b. (31) 

The ensuing fourth order differential equation will then involve function coefficients whose 
singularity is no more singular than ■\. Specifically, let A(x) = x'^ix)^ where T(x) = 
05(35^ + 3bx + X^) 7^ 0. Upon multiplying both sides of the translated fourth order equation 
{O^S{x) — 0) by T(x) , one obtains an expression of the form 

^Si'\x) + ^^S^^\x) + ^S^\x) + ^S^^\x) + ^S{x) = 0, (32) 

A AAA 

where all the C„(x) arc polynomials in x- 

Let us now multiply by x^i where p > 0. Upon doing so, through the operational 
procedure corresponding to 'integration by parts', we can rewrite the expression 

x'T{x)O^S{x) = ( E U.ix'^^ x"-'; p))S{x) + a,(T(x; S{x), ■ ■ ■ , S^'^x))) = 0, (33) 

j=0 

involving the total derivative expression noted, and functions of the power expressions, x'^^ 
where the exponent may be negative. 

Since the physical solution must decay exponentially, along the asymptotic directions 
[x ioo), wc can now integrate along a contour that sits on top of the real axis, except 
near the origin where it deviates around it. Denote this contour by C. Define the moments 
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(34) 



If p > 0, because of the analyticity of the solution, these moments correspond to the Ham- 
burger moments taken along the real axis: /ip = dxx^'^ixjj for p > 0. 
The moments will satisfy the moment equation (i.e. p — p — 2 > —2) 



3Es, 



+Pp-i 



p(-4 - p + 4p^ + p^) 



2Es7 



rGEso 



■p{p + 2) (6 + 12b^C2E + 7p + p^- m^s) 



62 



-2c2^(4 + 5p + p^) + 556^(14 + 25p + 8p^)) 



I2ES2S5 



63 



{-2c-2E{l2 + 8p + p^) + 6^55(168 + 169p + 38p^)) 



(42 + 30p + 5p2) +;ip^3 



+Pp+4 



2Esi 



63 



18S^s^ + 6s5(56 + 31p + 4p^)) 



64 



(162^^5^ + 655(42 + 19p + 2p')) + 



432^3^: 



3^3. 



144^3, 



67 



36^3< 



65 



r324£;3s^ 



68 



4^3 <,3 



66 



69 



0, (35) 



where p > —2. 

We note that for p > 0, the first nine Hamburger moments, {po, • • • , A^s} generate all of 
the other moments {/ip, p > 9) through the linear relation 



(36) 



where the coefficients Mp^^ satisfy the moment equation (Eq.(35), for p > 0) with respect 
to the p-index. In addition, they satisfy the initialization conditions M^-^^g^ = 81^^^^, for 
< ^1,2 < 8. 

If the moment equation is restricted to p > 0, then the associated differential equation 
would be of the form Oj^S{x) = ^(x); with V an unknown distribution supported at the 
origin (as discussed before). By also including p = —1, —2, we are imposing V — 
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For these two values of the p index, the moment equation becomes {p — —1) 



1^0 
+1^3 



y2Es2 
2Es2 



[-10c2E + 37b^S5] 



204ES2S5 



+ A*2 



162E^S2 + 25bs5) + fa 















[ ^5 J 




[ ^6 J 


+ 1^6 



ISE^si + 296S5) 



36^3^^ 



&8 



4Eh^, 



67 



E(//_i,//_2,/[/-4), 



(37) 



where 



Es \8Es 
E(/x_i,//_2,/X-4) = (l8bEs2S5iJ,_i j^{12b'^C2E - I2b^ 85)11-2 j-^jJi-A), 



and [p ~ —2) 

-24^S2S5 



+ ^1 



324^3^^ 



6^59 



66 



18£;^s^ + IO6S5) +//2 

3 
2 



(38) 



144^^5? 



+ A*6 



r2£;s, 



64 

36£;3si 



162£;^s| + I26S5) +//3 



432£;3s3 



68 



+ /"7 



4^'^s 



3^3. 



65 



69 



-E(yU_i,/X_2, At-4)- (39) 



Therefore, we can use both relations to constrain the /ig moment in terms of the moments 
{//o, A*i, ■ ■ ■ , A*?}- These in turn can be used with the moment equation (Eq.(35)), for p > 0, 
to generate all of the other Hamburger moments. Thus, we have 



defined by subtracting Eq.(39) from Eq.(37) (after dividing the former by |). 

Incorporating this within the expression in Eq.(36), 

7 _ _ 7 _ 



(40) 



(41) 



yields 



where 



M 



i=0 



Ms,e, p = 8 



(42) 



M, 



p,£ 



(43) 



and Mf^^^j = ^^i/j; < ^1,^2 < 7. This is the moment generating relation required in 
implementing EMM. 
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An additional important point is specifying the choice of normahzation. Since we are 
only dealing with Hamburger moments, some of the odd order moments may be negative. 
However, the even order moments are positive, and by working with the first 8 of them 

ue = i^2e, (44) 

we can impose the normalization 

E«^ = l. (45) 

e=o 

Working with the positive quantities, Ui, makes the implementation of the linear program- 
ming algorithm more efficient (since one is working within a bounded polytope corresponding 
to the unit hypercube, [0, 1]^). 
Thus, we take 

7 

Ue = IJ,2£ = Y1 ^2^A/"^. ' (46) 
for < i <7. The inverse of this relation is denoted by 

7 

/^^= E (47) 
ev=o 

and substituted in Eq.(42) producing 

7 

f^p^ Mp,eNe,e^Ui^, (48) 
e,ev=o 

or 

7 

I^P = YMp,eUe, (49) 



for p>0, where Mp,e = EL=o ^p/^,^^^„^ 

Finally, the normalization condition is incorporated {uq = 1 — J2l=i ui), producing 



7 

Atp = E -^P/^^' (50) 



where 
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Mo,e, p = 



Mp,e - Mo,e, P>1 



(51) 



and 



1,£ = 



(52) 



B. Numerical Results 



The numerical implementation of the EMM procedure is given in Table I. We only 
investigate the lowest energy state. The parameter Pmax denotes the maximum moment 
order generated, {/Xp|0 < p < Pmax}- We see that as the rotation angle increases, the 
tightness of the bounds decreases, as expected. Furthermore, for ^ > ^T'ao EMM solution is 
generated, consistent with the fact that the bound state becomes exponentially unbounded, 
in the corresponding direction in the complex x plane. 
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TABLES 

TABLE I. Bounds for the Ground State Energy of the —ix^ Potential Using Eq.( ) {Pmax = 40) 



El;o Eu-o 

0* 1.15619 1.15645 

.01 1.15617 1.15652 

.05 1.15510 1.15761 

.10 1.14201 1.17298 
.15 1.0 1.3 
.20 .5 8.0 



*Using Stieltjes EMM formulation by Handy(2001) 
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